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Abstract 



The spectrum and the mixing angle of the light scalar nonet mesons are stud- 
ied using the extended Nambu-Jona-Lasinio (NJL) model as well as the improved 
ladder approximation (ILA) of QCD with C/a(1) symmetry breaking interaction. 
The ILA is an approximation that is consistent with chiral symmetry and con- 
sists of the rainbow approximation of the Schwinger- Dyson equation and the 
ladder approximation of the Bethe-Salpeter equation. Improvement is made by 
the use of the running coupling constant according to the Higashijima-Miransky 
method. The C/a(1) breaking is supposed to come from the coupling of light 
quarks to instanton, and is represented by a contact six-quark vertex, called the 
Kobayashi-Maskawa-'t Hooft (KMT) interaction. The strength of the KMT in- 
teraction in the NJL model is determined so as to reproduce the electromagnetic 
decays of the rj meson. That in the ILA approach is determined so as to re- 
produce the pseudoscalar meson spectrum. This interaction is so strong that it 
causes large mixing of (uu + dd) and ss, in the 77 meson. 

In the extended NJL model, we study the qualitative features of the scalar 
meson spectrum. In the scalar nonet spectrum, the KMT interaction is found to 
give the right ordering of a — a masses and a few hundred MeV mass difference 
between the a and ao mesons. We also find that the strangeness content in the 
a meson is about 15%. 

In the ILA approach, we confirm the qualitative features of the results from 
the extended NJL. We obtain the mass spectra of the light pseudoscalar nonet 
mesons and the a, ao mesons which are consistent with current experimental data. 
They are both reproduced with the same parameter choice. On the other hand, 
we show that the obtained f and state can not be identified directly with 
the experimental data. It may suggest that the / (980) and k(700 — 900) states 
may not be explained only by the quark-anti quark state. 

We also find that the strangeness content in the a meson is about 5%. We 
obtain the result that the Ua(X) breaking interaction reproducing the mass spec- 
trum of the pseudoscalar meson gives large effects on the scalar meson mass 
spectrum. 
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Chapter 1 
Introduction 



Understanding low-lying hadron spectrum is one of the most challenging problems 
in the quantum chromodynamics (QCD). The spectrum is highly nontrivial due 
to the nonperturbative complexity, such as dynamical chiral symmetry breaking 
(D^SB), axial U(l) anomaly. For instance, the pseudoscalar mesons are off-scale 
light, if we suppose that they are bound states of a quark and an antiquark, while 
their spin-flip partners, i.e., vector mesons, are normal with masses about 2/3 of 
the baryon masses. This "anomaly" in the pseudoscalar mesons is attributed to 
their Nambu-Goldstone-boson nature associated with the D^SB. This is strongly 
in contrast with the heavy meson spectroscopy, such as heavy quarkonia and 
heavy flavor mesons. There the spectrum is more like the hydrogen atom with 
slightly stronger fine and hyperfine splittings. 

It should be noticed that the low-lying hadrons are the key to explore the 
complicated QCD vacuum, as in QCD we are not able to "measure" the bulk 
properties of the ground state, which can be accessed directly in the case of 
condensed matter physics. Thus it is important to explore the properties of the 
low-lying hadrons from the viewpoints of QCD dynamics and symmetries. 

Another nontrivial effect comes from Ua(1) symmetry, which is expected to 
be broken by anomaly. Weinberg showed that the mass of rj' should be less than 
V^m^ if C/a(1) symmetry were not explicitly broken. [1] Thus the C/a(1) symmetry 
must be broken. Later, 't Hooft pointed out the relation between C/a(1) anomaly 
and topological gluon configurations of QCD and showed that the interaction 
of light quarks and instantons breaks the Ua(1) symmetry. [2] He also showed 
that such interaction can be represented by a local 2Nf quark vertex, which is 
antisymmetric under flavor exchanges, in the dilute instanton gas approximation. 
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The dynamics of instantons in the multi-instanton vacuum has been studied by 
many authors, either in the models or in the lattice QCD approach, and the 
widely accepted picture is that the QCD vacuum consists of small instantons 
of the size about 1/3 fm with the density of 1 instanton (or anti-instanton) per 
fm 4 .[3] 

According to such instanton vacuum picture, the hadron spectrum shows its 
signature. The 77 — 77' mass difference is the obvious one, which can be understood 
by flavor mixing in the I = ^ (uu + dctj and ss. Without the flavor mixing, 

(uu + dctj and ss would form mass eigenstates, and thus the ideal mixing is 
achieved. This is natural if the Okubo-Zweig-Iizuka (OZ1) rule applies. However, 
the OZI rule is known to be significantly broken in the pseudoscalar mesons. For 
instance in the Nambu-Jona-Lasinio(NJL) model, the electromagnetic 77 decay 
processes suggest that the mixing of ^ (uu + dctj and ss is indeed strong so that 
the 77 meson is close to the pure octet state. [4] 

Recently, the scalar mesons, J 71 " = + , attracts a lot of attention by two 
reasons. [5] (1) Experimental evidence for a (I — 0) scalar meson of mass around 
500-800 MeV is overwhelming. [6, 7] Especially the decays of heavy mesons show 
clear peaks in the nn invariant mass spectrum. (2) The roles of the scalar mesons 
in chiral symmetry have been stressed in the context of high temperature and/or 
density hadronic matter. [8] It is believed that chiral symmetry will be restored in 
the QCD ground state at high temperature (and/or baryon density). Above the 
critical temperature of order 150 MeV the world is nearly chiral symmetric and we 
expect that hadrons belong to an irreducible representation of chiral symmetry, 
if we neglect small mixing due to finite quark mass. The pion is not any more 
a Nambu-Goldstone(NG) boson, and has a finite mass and should be degenerate 
with a scalar meson, i.e., a. 

For Nf = 3, we expect to have a U(3) nonet of light scalar mesons, which are 
the chiral partners of the pseudoscalar nonet. We assume that a is the member 
of this nonet. The naive assignment of the flavor component of a (I = 0) and 
a (I = 1) is uu + dd and uu — dd, respectively. The masses of the uu + dd and 
uu — dd mesons must become degenerate, if the flavor mixing interaction is absent. 
However, the mass of ao(980), which is the lightest ao meson, is larger than that of 
a by about 400 [MeV]. We consider that this mass splitting is explained again by 
the Ua(X) breaking interaction. As is in the case of the light pseudoscalar mesons, 
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the Ua(1) breaking interaction contributes significantly for the light scalar meson 
nonet and then a — mass splitting is explained by the flavor mixing. The goal 
of this paper is to study the role of the £/a(1) breaking interaction to the mass 
spectrum of the light nonet scalar meson. Here we study the lowest scalar meson 
states in the / = 0, 1, 1/2 channels and the second lowest one in I = channel 
which we call a, a , Kq and f , respectively. The identification of these states 
with the experimentally observed states will be given later. 

We first employ a simple model, the extended Nambu-Jona-Lasinio(NJL) 
model. It is the simplest possible quark model with the correct symmetry struc- 
ture for the present purposes. [11] The chiral symmetry is broken both explicitly 
by the quark mass term and dynamically by quark loops, while the Ua(X) sym- 
metry is broken by the Kobayashi-Maskawa-'t Hooft (KMT) interaction. [2, 12] 
The KMT interaction is a contact six-quark interaction, which represents the 
Ua{1) breaking effects induced by the instanton configuration of the gluon in the 
vacuum. 

After understanding the qualitative features of the scalar meson spectrum in 
the extended NJL model, we proceed to the analysis using the improved lad- 
der approximation (ILA) [14] of QCD with the KMT interaction. The ILA of 
QCD is the approximation in which the gluon is integrated out and is written 
only in terms of the quark degrees of freedom. The ILA is an approximation 
that is consistent with chiral symmetry. We employ the rainbow approximation 
of the Schwinger-Dyson (SD) equation for the quark propagator and the lad- 
der approximation of the Bethe-Salpeter (BS) equation for the quark-anti quark 
bound state. It is important to note that these two approximated equations are 
consistent with the chiral symmetry. Improvement is made by the use of the 
running coupling constant according to the Higashijima-Miransky [13] method. 
The KMT interaction is introduced in the ILA as an explicit interaction kernel. 
Naito et al. [17, 18] showed explicitly that the chiral symmetry is indeed pre- 
served in this formulation and properties of the pseudoscalar meson, such as the 
masses and decay constants, are consistent with those of the Nambu-Goldstone 
(NG) bosons of chiral symmetry breaking. They furthermore considered effects 
of the Ua(X) breaking and explicit chiral symmetry breaking due to the finite 
quark masses. They successfully reproduced the realistic mass spectrum of the 
pseudoscalar mesons. 
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In this thesis, we consider the scalar mesons in the same formulation. The 
advantage of this approach is consistency with chiral symmetry. It is found 
that the properties of the light pseudoscalar mesons are reproduced according 
to the expectations from chiral symmetry. We then treat the scalar mesons and 
the pseudoscalar mesons as the chiral partners. In this study, at first we fix the 
parameters so that the mass spectrum and the decay constants of the pseudoscalar 
mesons are reproduced and then we apply those for the scalar mesons. 

This subject has been studied in other approaches by many authors. The 
most direct approach might be the lattice QCD simulation. At present, however, 
it is difficult to compute the light singlet scalar meson. There have been studies 
by using the NJL model [19, 21] or the linear sigma model [22]. In those studies, 
the effects of the Ua(1) breaking interaction are weak so that the flavor mixing 
in the a meson is negligible. However, recently Takizawa et al. [4] showed 
that the Ua{^) breaking effects must be much stronger in order to explain the 
electromagnetic decays of the rj meson. The main aim of this study is to explain 
the anomalous ordering of scalar meson can be explained by the strong Ua{1) 
breaking effects. Indeed we will show that the a — ao splitting and significant 
flavor mixing is induced by the strong C/a(1) breaking interaction. We will see a 
large mixing of ss component in the a meson in our approach. 

Other pictures of the light scalar mesons proposed so far include the multi- 
quark state or the meson molecule state [9]. There may exist the mixing of 
quark-anti quark state with these state. However, it is difficult to believe that 
these multi-quark states are the main component of the light scalar mesons. We 
will show that the scalar meson nonet can be explained purely as the quark - anti 
quark state with the {7a(1) breaking effects. 

In chapter 2, we explain the theoretical background for this study. In chapter 
3, we present the extended NJL model analysis. In chapter 4, we present the 
formulation of the ILA approach. In chapter 5, we show our results and give 
discussions on the mass spectrum as well as the mixings. In chapter 6, conclusions 
are given. 



Chapter 2 

Theoretical Background 

2.1 Quantum Chromo dynamics 

Quantum Chromodynamics (QCD) is the SU (3) nonabelian gauge field theory, 
which describes the strong interactions of colored quarks and gluons. A quark 
comes in 3 colors and gluons come in eight colors. Hadrons are color singlet 
combination of quarks, anti-quarks and gluons. 

The Lagrangian of QCD with the gauge fixing terms is 

+d^c a + g fabeA^c b )c e (2.1) 
F; v = d,Al - d v Al - g,f ahc A\Al (2.2) 

D, = d, + tg ^-A* (2.3) 

where go is the QCD coupling constant, and the f a bc is the structure constant of 
the SU(3) algebra. The tpj represent the quark fields in which the j = 1,2,3 is 
the color label. The A^ represent the gluon fields with a = 1, 8. The quantities 
c a are called ghost fields which are introduced by in gauge fixing procedure. The 
ghost occurs only in loops, and never appears as asymptotic state. The £o is the 
gauge fixing parameter. The limit £ 00 defines the Landau gauge and we use 
this gauge in the improved ladder approximation later. 

The most characteristic property of QCD, which supports that the theory 
describing the dynamics of the hadron is QCD, is the asymptotic freedom. From 
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the experiments of the deep inelastic scattering of the lepton from the hadron, it 
was found that the quarks are almost free in high momentum transfer region. This 
property is called asymptotic freedom and it was shown that only the nonabelian 
gauge field theory exhibits it. Then the extra symmetry of the nonabelian gauge 
field theory is identified with color symmetry which is suggested by the quark 
model. This completes QCD. 

To see the asymptotic freedom nature of QCD we show the running coupling 
constant a s (q 2 ) = g 2 R (q 2 ) / An . It is obtained by the one-loop renormalization 
group calculation as 



where Nf is the number of quark flavors with mass less than \fq*. The a s (q 2 ) 
is small at large q 2 . Consequently in the high energy region, the perturbative 
analysis of QCD makes sense and it is applied to the deep inelastic scattering, 
e + e~ annihilation process, decays of heavy quarkonia and so on. 

In contrast, the a s (q 2 ) becomes large as q 2 decreases and then the perturbation 
approaches become poor approximations. In this region, the nonperturbative 
effects, especially the confinement and the dynamical chiral symmetry breaking, 
are important. The confinement means that only colorless states are physically 
realized and hence quarks cannot be observed in isolated state. The dynamical 
chiral symmetry breaking is explained in the next section. Because the complexity 
of low energy QCD, various approaches have been tried such as to solve QCD on 
the discretized space time numerically (lattice-gauge theory), construct effective 
theories and models (th potential quark model, the bag model, the Skyrme model, 
the NJL model), or to consider some limit as N c — > oo. 

2.2 The Chiral Symmetry Breaking 

If one supposes that all the masses of Nf quarks are zero, the QCD Lagrangian 
satisfies U fL (N f )xU fR (N f ) = U v (l)xU A (l)xSU fL (N f )xSU fR (N f ) symmetries, 
ie., the invariance under the transformations 



®s{q 2 ) 



12tt 



(2.4) 



(33 - 2N f ) ln(g 2 /A 2 ) 



+ ••• 





(a = 0,...,N 2 f -l). (2.6) 



(2.5) 



2.2. The Chiral Symmetry Breaking 



7 



Although the quark mass term breaks those symmetries explicitly, they remain 
approximate symmetries of QCD for light flavors i.e., u, d and s. 

Among these symmetries, some of them are broken in the vacuum state. The 
£7a(1) symmetry is broken explicitly by the anomaly and we explain it in the next 
section. 

We here concentrate on the SUfi{3) x SUfn(3) chiral symmetry, which is 
broken spontaneously in the vacuum. Then the symmetries of the QCD vacuum 
are £/y(l) x SU /(3). In this section, we explain the chiral symmetry and its 
breaking in the case of Nf = 2 for simplicity. At first, we show it using two 
simple models, the linear sigma model and the two-flavor Nambu-Jona-Lasinio 
(NJL) model, which have the SU (2)l x SU {2)r chiral symmetry. Then we explain 
the chiral symmetry of QCD. 

The Lagrangian of the linear sigma model is 

£ = ipi^ip + ^<9 M 7T(9 M 7r + ^-d^crd^cr — g$(cr — ir ■ 7T7 5 )-?/> 

+ ^(- 2 + - 2 )-^ 2 + - 2 ) 2 (2.7) 

where ip is the doublet nucleon = T (p n) and 7r is the triplet of the pion fields 
77 = T {it l 7r 2 7T 3 ). This Lagrangian is invariant under the SU /l(2) x SU /_r(2) 
chiral transeformation, 

^l,r^^l,r = U l ,r^l,r X=(<r + iT-w)^V = U L EU} t (2.8) 
f4, R = exp(^y). (2.9) 

It should be noted that this symmetry is isomorphic to 0(4) symmetry, which is 
represented by the 4-dimensional (a, 7r) space. When u 2 is positive, the potential 
energy 

n-,-) = ^(- 2 + - 2 )-^ 2 + - 2 ) 2 (2.10) 

forms a 4-dimensional wine bottle shape. 

The classical minimum of V(a, 7r) gives the set of degenerate ground states 
satisfying 

a 2 + ir 2 = ^. (2.11) 
A 
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Let us choose a particular ground state, 

(a) = ^ = v, <tt) = (2.12) 

Then one sees that the chiral SU(2)fi x SU(2)fR symmetry is broken sponta- 
neously by the ground state (2.12). Defining b = a — v , we then rewrite the 
Lagrangian as 

C = i>{i$- gv)i> + l - [d^a - 2u. 2 b 2 } + \d^ ■ c^tt 

-g$(a - ir ■ 7r 75 )^ + \vb{a 2 + tt 2 ) - ^[(a 2 + tt 2 ) 2 - v% (2.13) 

This Lagrangian is not any more symmetric under the rotation in the (a, tt) space. 
One sees that the pion becomes massless, while the scalar meson b acquires a 
mass = y/2fjF, and also the fermion gets a mass M = gv. 

The appearance of the massless pion is a concrete example of the Goldstone 
theorem: if a theory has a continuous symmetry of the Lagrangian is sponta- 
neously broken, there must be a massless boson, which is called Nambu-Goldstone 
(NG) boson. In this case, the pion is the NG boson. The order parameter of the 
chiral symmetry breaking is v — (a). In this model, the pion decay constant f n 
defined by 



0\A^\n\p)) = tUp^ (2.14) 

coincides with v in the tree level. 

Next we show the two flavor and Nq color NJL model. This model is written 
in terms of quark degrees of freedom as 

C = WW + \g s ]T [(^rV) 2 + (^ 75 rV) 2 ] . (2.15) 
z i=0 

where ip = T {u d). This Lagrangian is invariant under the U(2)f L x U(2)f R 
chiral transformation 

1>l,r - i>' L , R = exp(^[ iR ^)^ i = 0,..,3. (2.16) 

In this case, the order parameter of the chiral symmetry breaking is the value of 
the composite operator (ipip) and then the symmetry is dynamically broken. 



2.2. The Chiral Symmetry Breaking 
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We analyze this model in the leading order of the 1/Nq expansion. Using the 
mean field approximation, the Lagrangian is rewritten as 



£mfa = ^9 + AW (2.17) 

(2.18) 

(2.19) 



M = - A (^) = | limtr[^(x)] (2.18) 



Eqs. (2.17) and (2.18) lead to the following equation for M: 

M-^/^j— (2.20) 
Rotating the integral into Euclidean space (p° — > ip A ), we find the equation: 



M = -nr I d P „9 , iT79 ( 2 - 21 ) 



4tt 2 7 p 2 + M 2 
where A is the ultraviolet cutoff. After the integral, we obtain 



M 3 A 2 + M 2 
M(k-1) = k— In— ±5- (2.22) 



4tt 2 



(2.23) 



There is always the trivial solution, M = 0, to this equation, and when k > 1, 
there is also a nontrivial one with M^O. In the case of the M^O, the value 
of is n °t zero an d it shows that the chiral C//l(2) x Ufn(2) is dynamically 

broken to the Uf v (2). The quarks acquire the dynamical mass M. 

According to the dynamical symmetry breaking, there should be the four NG 
bosons. To confirm this, we consider the Bethe-Salpeter (BS) equations for the 
pseudoscalar mesons and the scalar mesons. The BS equation in the present 
model in Euclid space is written as 

p P r d 4 k 

^F n ,m(Q + l^)Xnimi(Q, P)SF mim (q ~ ~) = J ^y^mn,n 2 m 2 {<i-, P)Xn 2 ,m 2 (Q', P) 

r J )nm(r 3 ) m2n2 + (i7 5 r J )„ m (i7 5 r J ) m2n2 (2.24) 



K =z^ 

1 ^-mn,n 2 m 2 « ^ 



where \ is the BS amplitude defined as 

Xnm(q;P) = e- iP ^ J d\x-y) (0\TMx)^ m (y)\P) e iq{x - y) ■ (2.25) 
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\P) is the bound state with momentum P. The BS amplitude for the pseudoscalar 
meson, x^ P \i an d the BS amplitude for the scalar meson, x^ S \i can De written 
in terms of four scalar amplitudes as 



xP n J<i;P) = ic T i 



~<P ( s P \q; P) + <P ( p\g; P)i + <P ( S\q; P)f+ 4 P) (?; p)\ 
4 s \q; p) + P)i + ^(?; p)?+ 4%; P)\m - if) 



where r a is the Pauli matrix which denotes the flavor structure of the meson 
state. In the chiral limit and in the leading order of the 1/Nq expansion, the BS 
equations becomes extremely simple. Furthermore we takes only the first order 
in P 2 /A 2 , then we obtain the following BS equation: 

(g 2 + ^)4 p V,n = ^- 2 fdee4 p \k 2 ,p 2 ) (2.26) 

(g 2 + Y + 4M )^^ p2 ) = 2dk2k ^s\k 2 ,P 2 ). (2.27) 

These equations are satisfied for each flavor j. 

If one identifies the amputated function x(q 2 ) = (M 2 — q 2 )x and substitute 
the ^ — > M 2 , Eq. (2.26) coincides with Eq. (2.20). Then we can easily solve it. 

The spectrum of Mp = —P 2 of pseudoscalar bound states and the spectrum 
of Mj = — P 2 of scalar bound states are determined from the following equations 

k-1 = k — — — In [ -5-2- 2.28 

M 2 -^, / A 2 + 3M 2 -^\ , x 

k-1 = « r^ln 2 - 29 

A 2 I 3 m 2 -^ J 



2 



We obtain for M ^ case 

M P = (2.30) 
M s = 2M. (2.31) 

The four pseudoscalar mesons are massless and they are the NG bosons for the 
chiral symmetry breaking £//l(2) x C//h(2) — > £//y(2). 

In the case of QCD, if all the u, d, s quarks are massless, which is called the 
chiral limit, the chiral SU /l(3) x SU /_r(3) symmetry is exact. Although it is 
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broken by the quark mass term explicitly, since those masses are much smaller 
than the typical hadron masses (~ 700[MeV/c 2 ] — l[GeV/c 2 ]), this symmetry is 
approximately satisfied. By considering the chiral symmetry and its dynamical 
breaking, many properties of light hadrons are understood. Especially, the fact 
that the pion is extremely lighter than the other mesons is explained as the pion 
is the Nambu-Goldstone (NG) boson accompanying the dynamical breaking of 
the chiral symmetry. Practically, in the QCD vacuum the chiral symmetry is 
dynamically broken and the residual symmetry is SUv{3)- In that time, eight 
symmetry is broken and then there are eight NG bosons. The light pseudoscalar 
octet mesons excluding n' are the NG boson. The rj and r( mesons are the mixed 
state of the octet NG boson state and the singlet state which is not the NG 
boson state. This non-zero value of the quark condensate can not be obtained 
by the perturbative calculation. It is caused by the nonperturbative effects. The 
ladder approximation which we use in this study is the convenient way to obtain 
the non-zero value of quark condensate. 

Additionally, in the chiral symmetry broken vacuum the effective quark mass 
is generated. 

2.3 Ua(X) Symmetry Breaking 

For understanding the properties of the light hadrons, it is necessary to consider 
the C/a(1) symmetry, which is expected to be broken by anomaly. The QCD 
Lagrangian is symmetric under the Ua(X) transformation. However, Weinberg 
showed that the mass of rj' should be less than y/Sm^ if C/a(1) symmetry were not 
explicitly broken. [1] The experimental value of the eta' mass is 957.78±0.14[MeV] 
and that of the pion mass is ~ 140 [Me V]. Thus the Ua(1) symmetry must 
be broken by the anomaly. Later, 't Hooft pointed out the relation between 
U A (l) anomaly and topological gluon configurations of QCD and showed that 
the interaction of light quarks and instantons breaks the Ua(X) symmetry. [2] 

The instanton is the classical topological solution of the gauge field in the 4- 
dimensinal Euclid space. Then it is self-dual F£ v = F* u where Fg v = \e^F^. 
The BPST [23] instanton solution in the singular gauge is given by 



(2.32) 
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e a »v p,v = 1,2,3 

Va^u = { Sap v = 4 (2.33) 
Sau P = 4 



where p is the instanton size. This solution has the non zero value 

q = I dAx ^ F ^ r (2 ' 34) 

This integral is the value related to the topology which called as Pontryagin index 
q. In the Minkowski space, the instanton gives the transition amplitude from the 
vacuum of topological index Q to the vacuum of topological index Q+q. 

The QCD vacuum may be considered as the multi-instanton state. The dy- 
namics of instantons in the multi-instanton vacuum has been studied by many 
authors, either in the models or in the lattice QCD approach, and the widely 
accepted picture is that the QCD vacuum consists of small instantons of the size 
about 1/3 fm with the density of 1 instanton (or anti-instanton) per fm 4 .[3] 

If the vacuum consists of instanton configurations, light quarks are affected 
by the instantons. Especially for a massless quark there exists a zero mode which 
satisfies 



-i)pip o = 0. (2.35) 



The solution is given by 



*,(«) = -t^^l^ * ( 2 .36) 



71 v^2 {x 2 + P 2 fl 2 

where is the hedgehog structure for the SU (2) color and spin. This solution 
satisfies the 75^0 = i>o and then the instanton induces transition from the left 
hand ipoi to the right hand ipoR m &y occur. Since the zero mode exists in each 
flavor, the Green's function the Nf quark and Nf antiquark (llN f '4 , R'4 , i)j, which 
is not Ua(1) invariant, does not vanish. 

This effect can be represented by an low energy effective interaction with 2Nf 
quark vertex. In the case of Nf = 3, it is given by a six-quark interaction as 

£ 6 = G D { det [^(1 - 75)^] + det [^(1 + 75)^] } (2.37) 

which is called the Kobayashi-Maskawa-'t Hooft (KMT) interaction. This in- 
teraction is not Ua(1) symmetric and is SU(3)fR x SU(3)fR symmetric. By 
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flavor 


pseudoscalar scalar 


octet 
singlet 


attractive repulsive 
repulsive attractive 



Table 2.1: The effects of the Kobayashi-Maskawa-t'Hooft interaction on the pseu- 
doscalar qq state and the scalar qq state. 

substituting the one ipip by the condensate or mass term , the four-quark flavor 
mixing interaction is produced. 

The effects of the KMT interaction on the scalar qq state and the pseudoscalar 
qq state are summarized in Table. 2.1. From this table, one sees in how the 
KMT interaction induces the r] — r)' mass difference, which can be understood 
by flavor mixing in the I = 4j (uu + dctj and ss. Without the flavor mixing, 
rj = (uu + ddj and rj' = ss would form mass eigenstates, and thus the ideal 
mixing is achieved. This is natural if the Okubo-Zweig-Iizuka (OZI) rule applies. 
However, the OZI rule is known to be significantly broken in the pseudoscalar 
mesons. By the flavor mixing effect of the KMT interaction, these states are 
mixed. Then the n approaches to the flavor octet state (uu + dd — 2ssj and rj' 
approaches to the flavor singlet state (uu + dd + ssj. For instance, according 
to the Nambu-Jona-Lasinio(NJL) model, the electromagnetic r] decay processes 
suggest that the mixing of (uu + ddj and ss is indeed strong so that the 
n meson is close to the pure octet state. [4] Then in 77, the KMT interaction 
works as the attractive force which competes with the effects of the increase 
of the ss component. On the contrary, in rj', the KMT interaction works as 
the repulsive force which competes with the effects of the decrease of the ss 
component. According to the improved ladder approximation approach which 
we will show later, as the KMT interaction increases, rj mass does not change 
after slightly increasing and n' mass increase monotonically. Consequently, the 
77 — 77' mass difference is explained by the Ua(1) breaking effects of the KMT 
interaction. 

In this study, we investigate the role of this C/a(1) breaking effect in the light 
scalar nonet meson mass spectrum. 
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2.4 The Improved Ladder Approximation of QCD 

The improved ladder approximation of QCD is based on the approximation in 
the Schwinger-Dyson (SD) and Bethe-Salpeter (BS) equation of QCD. The gluon 
degrees of freedom are integrated out and are represented by the quark-quark 
interaction kernels. Both SD equation and BS equation are restricted to the 
ladder diagrams of gluon exchange between quark lines and the gluon polarization 
and the vertex correction are taken into account through the running coupling 
constant approximately. 

The origin of this approach is the Higashijima-Miransky approach [13] for 
the Swinger-Dyson (SD) equation for the quark propagator, which keeps consis- 
tency with asymptotic freedom by using the running coupling constant of one 
loop renormalization group calculation. Since that running coupling constant 
diverges at low momentum, an infrared cutoff are introduced. Under the rain- 
bow approximation the SD equation in the Landau gauge shows the D^SB and 
gives the running quark mass consistent with the renormalization group at high 
momentum. 

Later, Aoki et al. [14] applied the same idea to the BS equation and shows 
that the combination of the rainbow SD and the ladder BS equation is consistent 
with chiral symmetry. It was shown that in the chiral limit, the pseudoscalar 
meson obtained as solutions of the BS equation have the expected properties 
consistent with chiral symmetry. 

The shortcomings of this approach may be summarized in the following four 
points. The first one is that the infrared cutoff can not be deduced from QCD. 
Although the strength of the D^SB depends on it strongly, it should be given 
phenomenologically. The second one is that results are gauge dependet. In the 
ILA, the Landau gauge is employed at the gluon propagator. To remove the 
gauge dependence, the higher order may be necessary. [15] The third one is that 
in ILA a large Aqcd is needed to reproduce the appropriate magnitude of the 
chiral symmetry breaking. It may indicate that the chiral symmetry breaking can 
not be explained only by the ladder approximation of the one gluon exchange. 
The last one is that the axial Ward identity of QCD is not satisfied by the non- 
local coupling constant of the gluon exchange interaction. However, by using the 
modified axial vector current the axial Ward identity can be satisfied. [16] In this 
study, this modified axial vector current is used. 
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A realistic model based on the ILA has been studied by Naito et al. [18] 
They introduced finite quark mass and further C/a(1) breaking interaction in the 
interaction kernel and calculated the pseudoscalar spectrum. 
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Chapter 3 

Extended Nambu-Jona-Lasinio 
model 



In order to clarify the roles of the breaking interaction in the scalar meson spec- 
trum, we first consider a simple model and analyze the pseudoscalar and scalar 
mesons. One of the simplest models which incorporate chiral symmetry and its 
dynamical breaking is the Nambu-Jona-Lasinio (NJL) model, which is extended 
by an additional Ua(1) breaking interaction. The Ua(1) breaking is considered 
to be caused by the interaction of quarks and instantons in the vacuum. Here 
the {/a(1) breaking interaction is represented by the six-quark interaction term 
called Kobayashi-Maskawa-t' Hooft (KMT) interaction. 



3.1 Formulation 

We work with the following NJL model lagrangian density extended to three- 
flavor case: 



C 

c 6 



Co + £4 + Cq, 

4> {id^ — rh) ip 
G f 



a=0 



G D { det [^(1 - 75)^] + det [^(1 + 75 )%] } 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



Here the quark field ip is a column vector in color, flavor and Dirac spaces and 
A" (a = . . . 8) is the Gell-Mann matrices for the flavor U(3). The free Dirac 
lagrangian £ incorporates the current quark mass matrix rh = diag(m u , ma, m s ) 
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which breaks the chiral Ul(3) x Ur(3) invariance explicitly. £4 is a QCD mo- 
tivated four-fermion interaction, which is chiral Ul(3) x Ur(3) invariant. The 
Kobayashi-Maskawa-'t Hooft determinant Cq represents the Ua{^) anomaly. It is 
a 3 x 3 determinant with respect to flavor with i,j = u, d, s. 

Quark condensates and constituent quark masses are self-consistently deter- 
minded by the gap equations in the mean field approximation, 

M u = m u -2G s (uu) -2G D (dd) (ss) , 
M d = m d - 2G S (dd) - 2G D (ss) (uu) , 
M s = m s -2G s (ss) -2G D (uu)(dd) , 



(3.5) 



with 



(qq) = -Tr^[iS q F (x = 0)] 
.Tr (c,D) 



A d 4 p 



(27T) 



p^ -M q + ie 



(3.6) 



Here the covariant cutoff A is introduced to regularize the divergent integral and 
means trace in color and Dirac spaces. 
The scalar channel quark-antiquark scattering amplitudes 

(P3,P4, out |pi,p 2 ;in) = (2tt) 4 5 4 (p 3 +p 4 -p x - p 2 )T qq - (3.7) 

are then calculated in the ladder approximation. We assume that m u = so 
that the isospin is exact. In the a and /o channel, the explicit expression is 



T- 
2 qq 



u(p 3 )\ 8 v(p 4 

u(p 3 )X°v(p 4 



A(q 2 ) B(<?) 
B(q 2 ) C(q 2 ) 



v(p 2 )\ 8 u(p 1 
v(p 2 )X°u(p 1 



with 



B(q 2 ) 



{2(G G 8 - G m G m )I°(q 2 ) - G 8 } , 



and 



detD(g 2 ) 
2 

detD(g 2 ) {~ 2 ( G ° Gs ~ G m G m )I m (q 2 ) - G m } , 
H -^{2(G G 8 -G m G m W)-Go}, 



G = ±G s -±(2(uu) + (ss))G D , 



G 8 
G m 



±G s -l((ss)-4(uu))G D , 



1 

371 



((ss) - (uu))G D 



(3.8) 

(3.9) 
(3.10) 
(3.11) 

(3.12) 
(3.13) 
(3.14) 



3.1. Formulation 



19 



The quark-antiquark bubble integrals are defined by 

- A d A p 



I°(q 2 ) = i 
/V) = i 
I m {q 2 ) = if 



(27T) 



r Tr^> D ) S F (p)X°S F (p + q)\° 



X r ( c -/' D ) 



T r ( c </> D ) 



S F (p)\ 8 S F (p + q)\t 
S F (p)\°S F (p + q)\* 



with q = pi + p 2 - The 2x2 matrix D is given by 

D(g 2 ) = 



D n (g 2 ) D 12 (q 2 ) 

£>2i (g 2 ) ^22 (g 2 ) 



with 



Dniq 2 ) = 2G I°(q 2 ) + 2G m I m (q 2 ) 

D 12 (q 2 ) = 2G I m (q 2 )+2G m I s (q 2 ) 

D 2l {q 2 ) = 2G 8 r\q 2 )+2G m I°(q 2 ) 

D 22 {q 2 ) = 2G 8 I 8 (q 2 ) + 2G m I m (q 2 ) 



(3.15) 
(3.16) 
(3.17) 

(3.18) 

(3.19) 
(3.20) 
(3.21) 
(3.22) 



From the pole positions of the scattering amplitude Eq. (3.8), the a-meson mass 
m a and the / - meson mass mj are determined. 

The scattering amplitude Eq. (3.8) can be diagonalized by rotation in the 
flavor space 



T- = 



xT, 



x 



u{p 3 )X 8 v{p 4 ) \ 
«(Ps)A°t;(p4) J 
( v(p 2 )X 8 u( Pl ) 

u(p 3 )\ a v(p 4 ) \ 
u(p 3 )\f°v(p 4 ) J 

^(p2)A <T u(pi) 



A{q 2 ) B{q 2 ) 
B(q 2 ) C(q 2 ) 



T -i 
J- a 



(3.23) 



D CT (g 2 ) 
D^{q 2 ) 



v(p 2 )X fo u( Pl ) ) ' 
with A 17 = cosflA 8 - sin#A°, A /o = sin#A 8 + cosflA and 

/ cos — sin 9 \ 
9 ~ sm9 cos9 



The rotation angle 9 is determined by 



tan 29 



2B(q 2 ) 



C(q 2 ) - A(qi) 



(3.24) 



(3.25) 



(3.26) 
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Note that therefore depends on q 2 . At q 2 = m 2 , 8 represents the mixing angle 
of the A 8 and A components in the cr-meson state. 

The Ua(1) breaking KMT 6-quark determinat interaction contributes to 
the scalar qq channel only by the form of the effective 4-quark interaction, which is 
derived from C§ by contracting a quark- ant iquark pair into the quark condensate. 
The explicit form of the effective KMT interaction is 



— ) (2 (uu) + (ss)) 



(^aV) 2 -(#W) : 



+ (ss) £ 



i=l 

7 



X 



- U'A^Ts 



i=4 



+ ( - ) (4 (uu) - (ss)) 



+ 



(^AV) -(^A 8 i 75 ^) 
^pj ((uu) - (ss)) [($AV) (^AV) + ($Afy) (V^AV) 



- (i>\°i^) + (V>A 8 n^) (^A°^75^)] I • (3.27) 



One can easily figure out from Eq. (3.27) that the Ua(1) breaking KMT inter- 
action gives the attractive force in the flavor singlet scalar qq channel. On the 
other hand, it gives the repulsive force in the isospin 1=1 (a ) and 7=1/2 
(Kq) channels. Because of the large strange quark mass, | (ss) | is bigger than 
| (uu) | , and therefore, the repulsion in the 1 = 1 channel is stronger than that in 
the 7 = 1/2 channel. 



3.2 Results 

We show our numerical results and give discussions on the mass spectrum as 
well as the mixings in this section. As the extended NJL model has been used 
in the analyses of the pseudoscalar mesons, here we have used the model pa- 
rameters fixed in the study of the electromagnetic decays of the rj meson. Since 
the rj meson properties depend on the strength of the Ua(1) breaking interaction 
rather sensitively, it is reasonable to determine the strength of the Ua(1) breaking 
interaction from the n meson properties. 
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The parameters of the NJL model are the current quark masses m u = m^, 
m s , the four-quark coupling constant Gg, the Ua(1) breaking KMT six-quark 
determinant coupling constant Gd and the covariant cutoff A. We take Gd as a 
free parameter and study scalar meson properties as functions of Gd- We use the 
light current quark masses m u = = 8.0 MeV to reproduce M u = ~ 330 
MeV (~ 1/3Mjv) which is the value commonly used in the constituent quark 
model. The other parameters, m s , G$ and A, are determined so as to reproduce 
the isospin averaged observed masses, m n = 138.0 MeV, mx = 495.7 MeV and 
the pion decay constant /„- = 92.4 MeV. When we take the different value of Gd, 
we go through the fitting procedure each time. 

We obtain m s = 193 MeV, A = 783 MeV, M u4 = 325 MeV, M s = 529 MeV 
and fx = 97 MeV, which are almost independent of Gd. The quark condensates 
are also independent of Gd and our results are (uu)^ = —216 MeV and (ss)^ = 
—226 MeV whenever we have fixed other model parameters from the observed 
values of m^, itlk and f n . 

We define dimensionless parameters, 

G c g = -G D (A/27r) 4 AV c 2 

Gf ee Gs(A/2tt) 2 N c . (3.28) 

As reported in Ref. [4], the experimental value of the 77 — > 77 decay amplitude 
is reproduced at about Gff = 0.7. The calculated r^-meson mass at G e § = 0.7 
is m ri = 510 MeV which is 7% smaller than the observed mass. G c § = 0.7 
corresponds to Gd{ss) /Gs = 0.44, suggesting that the contribution from £ 6 to 
the dynamical mass of the up and down quarks is 44% of that from £ 4 . The 
calculated value of T(r] — > 7r°77) is 0.92 eV at Gff = 0.7, which is in good 
agreement with the experimental data: T(r] — > 7r°77) = 0.93 ± 0.19 eV. 

Before going to present the numerical results for the scalar mesons, let us 
summarize the properties of the scalar mesons in the NJL model. In the SU L (2) x 
SUr{2) version of the NJL with no explicit symmetry breaking term, the cr-meson 
mass can be calculated analytically in the mean field + ladder approximation, i.e., 
m a = 2M U . The a meson is therefore regarded as the lowest bosonic excitation, 
whose mass is twice of the gap energy, associated with chiral symmetry breaking. 
It should be noticed that there is a cut above q 2 = AM 2 in the complex g 2 -plane 
of the quark-antiquark scattering T-matrix, which corresponds to the unphysical 
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decay: a — > qq. This is one of the known shortcomings of the NJL model. If one 
introduces a small symmetry breaking term, i.e., the current quark mass term, 
rricr moves up and gets the imaginary part corresponding to the a — > qq decay. 
[25] The pole position is in the second Riemann-sheet of the complex g 2 -plane, 
as is the case of ordinary resonances. It means that the Argand diagram for the 
T- matrix makes a circular resonance shape in the scalar qq channel. 1 

It should be noted that the physical decay mode of a, i.e., a — > tttt is neither 
taken into account in the ladder approximation. As this decay makes the a width 
significantly large, our result for o mass is qualitative rather than quantitative. 
Nevertheless, the results shown below show that the scalar mesons in the NJL 
model is realized as the chiral partner of the Nambu-Goldstone bosons, and that 
they give systematic behavior for the orders of the masses and the splittings. 

Let us now discuss our results of the scalar mesons. The calculated results of 
the scalar-meson masses, qq decay widths and the mixing angle 9 are shown in Fig. 
3.1, Fig. 3.2 and Fig. 3.3, respectively. The qq decay widths of the scalar mesons 
shown there are unphysical ones. We present them just for showing the pole 
positions in the complex g 2 -plane. When Gff is zero, our lagrangian does not 
cause the flavor mixing and therefore the ideal mixing is achieved. The a is purely 
uu+dd, which corresponds to 9 = —54.7°, and is degenerate to the ao in this limit. 
When one increases the strength of the {7a (1) breaking KMT interaction, the qq 
attraction in a increases and a state moves from the ideal mixing state toward the 
flavor singlet state. It means that the strange quark component of a increases as 
becomes bigger. Since the increase of the attractive force compensates with 
the increase of the strange quark component, m a is almost independent of the 
strength of the Ua(1) breaking interaction. The qq decay width of the a meson 
is very small, i.e., less than 2 MeV and therefore we neglect it in our calculation 
of the mixing angle. At = 0.7, the calculated mixing angle is 9 = -77.3°, 
corresponding to about 15% mixing of the strangeness component in a. 

Hatsuda and Kunihiro have discussed the masses and mixing angle of the 
isoscalar nonstrange (o^s) and strange (as) scalar mesons using the similar 
model. [19] They have reported a rather small mixing between a^s and as- The 
reason of the difference between their result and our result is the strength of 

1 The situation is quite different in the case of the vector meson. In the nonrelativistic limit, 
the scalar meson channel corresponds to the p-wave quark-antiquark state whereas the vector 
meson channel corresponds to the s-wave quark-antiquark state. See Ref. [26]. 
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Figure 3.1: The calculated scalar meson masses as functions of the effective 
coupling constant G e ^ of the C/a(1) breaking KMT interaction. The solid, dashed, 
dotted and dash-dotted lines represent m a , m ao , rriK* and to/ , respectively. 

the U A {1) breaking KMT interaction. The strength of the U A {1) breaking KMT 
interaction used in the present study is much stronger than that used in their 
study. They have determined the strength from the r( mass, while we have fixed 
it from the radiative decays of rj. Strong Ua(1) breaking interaction suggests 
that the instanton liquid picture of the QCD vacuum. [3] In Ref. [19], they have 
discussed the origin of the difference of the mixing properties of the scalar mesons 
and pseudoscalar mesons. We agree with their qualitative discussion, namely, the 
flavor mixing between a and f is weaker than that between rj and r] f . Shakin 
has also pointed out that the KMT interaction mixes the a^s an d cr 5 , while he 
assigned the lowest I = qq state to / (980). [20] 

Let us turn to the discussion of the a and Kq mesons. As shown in Fig. 3.1, 
both m aQ and m K * increase as increases. The slope for m ao is steeper than 
that for rriK*, which is consistent with the simple argument based on the form 
of the effective interaction Eq. (3.27). At G*§ = 0.7, the calculated masses are 
m ao = 816 MeV and m K * = 1002 MeV, therefore the C/a(1) breaking interaction 
pushes up the a and Kq masses about 161 MeV and 88 MeV, respectively. 
Although the effect of the Ua(1) breaking interaction on the Kq meson is smaller 
than that on the a meson, our numerical results show that it is not enough to 
support the existence of the light k state. 
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Figure 3.2: The calculated qq decay widths of the scalar mesons as functions of 
the effective coupling constant G c § of the Ua(1) breaking KMT interaction. The 
solid, dashed, dotted and dash-dotted lines represent m a , m ao , m K * and m/ , 
respectively. 

As for the /o meson, we have shown our results in Figs. 3.1 and 3.2. At 
G c p = 0, the /o state is expected to be pure ss state in our model. Because of 
the qq decay width, we cannot calculate the mixing angle for / . The calculated 
mass of the /o meson at G*§ = is m/ = 1.163 GeV which is above the ss 
threshold 2M S = 1.113 GeV. As shown in Ref. [21], the symmetry breaking effect 
by the current quark mass term pushes up the scalar meson mass above the 
qq threshold and the following relation is obtained by using the bosonization 
technique with the lowest order derivative expansion in the NJL model. 

( m scaiarmcson ~ (<?<? threshold energy) 2 ) oc m currcnt quark (3.29) 

Our results at G c £ = are m 2 - 4M 2 = 0.008 GeV 2 , m\, - (M u + M s ) 2 = 
0.060 GeV 2 and m 2 o - 4M 2 = 0.115 GeV 2 , respectively. The above simple mass 
relation therefore holds in our case too. Fig. 3.1 shows that the f meson mass 
is almost independent of the strength of the KMT interaction. The situation is 
just opposite to the a case, i.e., the increase of the repulsive force by the KMT 
interaction compensates with the decrease of the strange quark component of the 
f meson when one increases the strength of the KMT interaction. 

It should be noted here that in the SUl(3) x SUr(3) version of linear sigma 
model, not only the three-meson flavor determinant term but also the chiral 
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Figure 3.3: The calculated mixing angle of the a meson as a function of the 
effective coupling constant G c § of the Ua(1) breaking KMT interaction. 



invariant four-meson terms give rise to the a — a mass difference. [22, 9] We 
note that the extended NJL model does not give such type of interaction. 
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Chapter 4 



Formulation of the Improved 
Ladder Approximation Approach 



In this chapter, we present the formulation of the Improved Ladder Approxima- 
tion (ILA) of QCD. The approximation is applied to the Schwinger-Dyson (SD) 
equation for the quark propagator and to the Bethe-Salpeter (BS) equation for 
the pseudoscalar and scalar meson. 

4.1 Lagrangian 

The ladder approximation is the lowest order in the perturbation theory for the 
interaction kernels employed in this approximation may not be valid except for a 
heavy quark systems, where the coupling constant is small for this approximation. 
It, however, can be improved at low mass region by including a certain set of 
higher-order diagrams which bring the running coupling constant, a,s{p 2 )- 

In the present study we employ the following Lagrangian density the ILA of 



In (4.1), m denotes the bare quark mass matrix m = diag(m g , m g , m s ) where 
we assume the isospin invariance. £ge denotes the gluon exchange interaction 



QCD, 



ip = (u, d, s) T . 



(4.1) 
(4.2) 



2 Jpp', 




x^ m b)^(p , )^(?)^(9 , )e- j(p+p ' +<?+9> (4.3) 
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where 



/C™»''™W, q, q') = i?DV (^f - q -^j ( 7 ,T f 



) mm \ lu T a y 



(4.4) 



To reduce the expressions we use an abbreviation f (27r) 4 through this chapter. 
The indeces m,n, ■ ■ ■ represent the components of the color, flavor and Dirac 
space. In the gluon propagator we employ the Landau gauge, 



/ k-f M k v \ — 1 



(4.5) 



and the Higashijima-Miransky type running coupling constant g 2 defined as fol- 
lows. 



9 (Pe, 4) = 0(Pe ~ Qe)APe) + 0(q E - Pe)AQe) 



(4.6) 



with 



Ape) = 



i i 

2Po (1+tlF) 2 

1 3tiF-<o+2 
I Wo (1+*if) 2 



Stip — to + 2 — 



(t-to) 2 

*IF— to 



A 2 



A) 



1 llN c -2N f 



(4tt) 2 3 

Here, pe and denote the Euclidian momenta defined by 



for t < t < tiF (4.7) 
for t < t 

(4.8) 
(4.9) 



P = (Po,p) 

2 2-2 

P =Po~P 



Pe = (ip4,p) 

2 2-2,2 
P E = ~P =P +Po 



(4.10) 
(4.11) 



In eq.(4.7) the infrared cut-off tip is introduced. Above tip 9 2 {p 2 e) develops accord- 
ing to the one loop solution of the QCD renormalization group equation, while 
below t , g 2 {p 2 E) is kept constant. These two regions are connected by a quadratic 
polynomial so that g 2 {p 2 E ) becomes a smooth function. Here N c and Nf are the 
number of colors and active flavors respectively. We use N c = Nf = 3 in this 
study. The behavior of the running coupling constant is shown in the Fig. 4.1 to- 
gether with the one-loop QCD running coupling constant for A QCD = 600 [MeV]. 
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Figure 4.1: The running coupling constant of the QCD and the ILA 



£kmt is the Kobayashi-Maskawa-'t Hooft (KMT) interaction given by 



'KMT 



g Uflt t W \U Xl , U yi , U X2 , Uy 2 , U X:j , Uy 3 ) 

x { (xi)^h (yi)\ (y2)} MiPk (ys)] 

+3 [i) 91 {x 1 )tl) h (yx)} [i> g2 (x 2 )75^ /2 (y 2 )} bP 93 (a: 3 )75^/ 3 (ys)} } 



(4.12) 



where f±, g±, • • ■ are flavor indices, e denotes the antisymmetric tensor with e uds = 
1 and the symbol * at the end of the equation means to take the limit x±,X2, - ■ ■ — > 
x after all derivatives are operated. This interaction causes the Ua(X) symmetry 
breaking and also flavor mixing. 

We introduce a weight function w which is necessary so that the KMT inter- 
action is turned off at the high energy region. Then the asymptotic behavior of 
the ILA are kept. We use the following Gaussian form 



/02 + Q2 +d 2 + Q2 +d 2 +d 2 \ 

w(d Xl ;d yi ;d X2 ;d m ,d X3 ; dj = w " + " + * 3 . yi + y2+ VA 4.13) 
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w(fi 2 ) = exp(-K/x 2 ) (4.14) 

This weight function is convenient for numerical calculations as it is factorized as 

w(-p 2 - q 2 ) = w(-p 2 )w(-q 2 ). (4.15) 

The parameter k is taken as 

k = 0.7[GeV~ 2 ]. (4.16) 

This value corresponds to the form factor of the instanton of the average size p, 
about l/3[fm]. The instanton form factor, 



X E + P 2 P" 

can be identified for small x\ with the Fourier transformation of the weight 
function 



ocl-% + 



(4.17) 



X 



F.T.^)=Cexpl-^]ocl-^ + 



with 



Ak 



Ak = p 2 . 



Ak 



(4.18) 



(4.19) 



g(q 2 )Wnraij<^(9 2 ) + u 




u d 

Figure 4.2: The interactions of the ILA. 

Finally we introduce the ultraviolet momentum cutoff Auv 5 which is necessary 
to regularize loop integrals in the ILA approach. It is useful to introduce the cutoff 
function 



1 



f(p 2 



1 P 2 <A UV 
p 2 >A uv 



(4.20) 



4.3. Effective Action 



31 



By adding it in the kinetic term as 

$f(d 2 )(i$-m )iJ (4.21) 

all integrals are regularized. We choose a very large A uv (~ 100[GeV]) so that 
the results are not sensitive to the value of Auv- 

It should be noted that because the cut off function, f and w, depend on the 
momentum, the Noether currents are modified from the bare QCD currents. 



4.2 Renormalization of Quark Masses 

In our approach, the parameters of the Lagrangian are renormalized at the cutoff 
scale Auv- Since the Auv should be taken enough large, it may not be appropriate 
to fix the quark mass at such high energy. Instead of the quark mass m qo , m s0 at 
the scale Auv, we choose renormalized the quark mass at the momentum scale 
/j, = 2[GeV] as m qR = 5[MeV] and m sR = 100 [Me V] and calculate the m qo and 
m s0 in the following way. 

The renormalization constants Z mq and Z ms are defined by 

m q = Z m q m qw m s = Z m] m sR- (4-22) 



We take the renormalization condition as 

dB q {^) 



dm qR 



= (4.23) 

m qR 



dm 



sR 



= 0. (4.24) 



We define the quark condensate as follows 



where 



<: iKO)^(O) :>= - / Tr[^ ull (p)] + / Tr[S h F ave ] (4.25) 

4i — Tflf! 

Tr[Sr c (p)] = N C / d f± ~ 2 (4.27) 
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4.3 Effective Action 



To derive the Schwinger-Dyson (SD) equation and the Bethe-Salpeter (BS) equa- 
tion, we use the Cornwall- Jackiw-Tomboulis (CJT) effective action formulation 
[29]. 

The CJT effective action is defined as 



T[S F ] = tTrLn[S F ] - iTr[S^S F ] + T loop [S F ] 



(4.28) 



The last term of Eq.(4.28) is the residual term, and is given by the sum of all 
the Feynman amplitudes of corresponding to the two particle irreducible vacuum 
diagrams with two ore more loops. In these amplitudes, quark propagator is 
substituted by the full propagator 



S F (x,y) = {0\Ti/>{x)$(y)\0). 



(4.29) 



In this study, we employ the lowest loop order and the leading 1/N C contribution 
for the r loop [SV]. Then the CJT action becomes 

F[S F ] = iTrLn[S F ] - iTr[S^S F ] + T GE [S F ] + T KMT [S F ] 



(4.30) 



where 



- l - J d 4 xJC m ^ n ™(id xl ,id X2 ;id yi ,id y2 ) 

x (SFm 2 mi(z2,Zi)S Fn2ni (y 2 ,Z/i) - S Fm2Tll (x 2 , yi)S Fn2mi (y 2 , i4$l) 



KMT 



[5, 



r r ( f) 2 4- f) 2 4- f) 2 



Q2 + Q2 + Q2 +d 2 + Q2 +d 2< 



x (-tr(° c ) [S Ffl9l ( Vl , x 1 )]tr(° c ) [S Ff292 (y 2 , x 2 )]tr^ [S FJaga (y 3 , x 3 )} 

-3tr( D °) [S F fl91 ( yi , Xl )}trW [ l5 S F hg2 (y 2 , x 2 )]tr^ [ l5 S F hgs (y 3 , x 3 )\ 

(4.32) 

It should be noted that the global SUl(3) x SUr(3) symmetry is preserved 
within this approximation if the quark masses m q and m s vanish. 

In fact, the total effective action T[S F ] is invariant under the infinitesimal 
global chiral transformation 

S F (x, y) -> (l + *7sy £ a ) Sf{x > y) f 1 + tl5 Y 6a ) ' (433) 
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1 -1 



SrV- = — > — + ^s?)4 




F> 



Figure 4.3: The SD equation 



+ 




except for the quark mass term. 

From this truncated effective action, the rainbow SD equation and the ladder 
BS equation are derived consistently. 

4.4 Schwinger Dyson equation 

The Schwinger- Dyson equation is derived by the stability condition of the CJT 
action 

ST[S F ] 



5S F (x,y) 

Introducing the regularized propagators 



0. (4.34) 



S?(q) = f(q 2 )S F (q) (4.35) 
S?(q) = f(q 2 )S (q), (4.36) 



the SD equation in the momentum space becomes 

iR~l{„\ •nil-l/ >, C F f 1 -2/2 2\ 



xtr( DC )[^(p) 9i/ Jtr( DC )[^(A;) 92/2 ] (4.37) 
with the coefficient from the color space, 

N c 2N C k ' 

This equation is shown diagrammatically in Fig. (4.3). 
Generally the quark propagator is parametrized as 
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In the Landau gauge, it can be shown that the solution satisfies A(—q 2 E ) = 1. 
Then the SD equation becomes the integral equation only of the mass function 
B(q 2 ) and reads 

B q (- q >) = m q + f V *W, *■) {,<*■ - f) + * W - *■)} 

+ 8tt* Jo ( } k* + B q (-k*) Jo d/W(/) /2 + 5s (-/ 2 ) 

(4.40) 



B s {-q 2 ) = m s 



(4.41) 

4.5 Bethe-Salpeter equation 

To treat the pseudoscalar mesons and the scalar meson as quark-antiquark bound 
state, we use the homogeneous Bethe-Salpeter equation. It is derived by 

xq , S ^! S q F 1 , , . Xn'MMPs) = (4.42) 
db Fmn {x,y)db Fm , n ,{y',x') 

where 

Xn>,m>{y',x'-P B ) = (0\Tijj n ,(y')if; m ,(x')\P) (4.43) 

denotes the BS amplitude. The normalization condition is taken as 

{P B \P' B ) = (2n) 3 2P B0 5 3 (P B - P' B ) (4.44) 

and P B = \J M B + P\ where P B is the on-shell momentum of the meson. 
Introducing the regularized BS amplitude by 

xS >m '(?; Pb) = /(-£)x»',m'(?; P B )f(-d) (4.45) 

the BS equation in momentum space becomes 

sf 1 ( q+ )j ( R (q;P)sf\q.) 



4.5. Bethe-Salpeter equation 
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X\-- + 




Figure 4.4: The BS equation 



-iC f 



g 2 (-q 2 , -k 2 )iD^{q - k) 7 , X R (k; P)% 



Ik f(-k 2 + )f(-k 2 _) 
Jp 

x^[S% h (p)} {(l 5 )tT^[ l5X f'Ak;P B )} + (l)tr^)[ X ^(A;;P B )]}4.47) 



P, k f(-p 2 )f(-k 2 + )f(-kl) W { ^ ^ ^ 2* 



which is shown diagrammatically in Fig. (4.4). 

For the pion, the BS amplitude can be written in terms of four scalar ampli- 
tudes as in Ref. [18], 



A" 



xtLfc P) = icy [Mk; P) + Mk; P)f + <j> Q (k; P)F+ Mb p )^ ~ W) ^- 48 ) 

where A a denotes the flavor structure of the pion state. The neutral pion, for 
instance, is given by a = 3 



1 
0-10 




(4.49) 



On the other hand, for the r] and rf mesons, the BS amplitudes are written in 
terms of eight scalar amplitudes, 



loT p ) + p W + p ) p + P)\ 



2 
X s 



75 



P) + <p s P (k; P)^ + <P s Q (k; P)f+ Mk; P)\(ft - ftfl^ 



where the flavor matrices X q and X s are defined by 












( 









X q = 


1 


!)■ 


X s = 










(4.51) 


\o 









v0 





V2) 





Because of the flavor mixing effects caused by the KMT interaction, the q and s 
components of the BS amplitude are mixed. We identify the ground state solution 



36 Chapter 4. Formulation of the Improved Ladder Approximation Approach 

with the i] meson state and the first excited state with the rj meson state. The 
explicit form of the BS equation is given in the appendix. 

Similarly the pseudoscalar mesons, the BS amplitude of the scalar meson is 
parametrized as 

xlJk- P) = icy (W; P) + Mk; P)ff + P)F+ Mk; P)\{ft- m)-52) 

for a and 

xL = icy (f s (k\ p) + p)t + m-, nr+ p)\m - m) 

+ i C y (^(fc; P) + Mk; P)t + <i> s Q (k; P)F+ Mk; P)\{ft ~ W^3) 

for o or /o, respectively. 

In the numerical computation, we treat the BS equation in the Euclidean 
momentum region. Then the physical solution, which corresponds to negative 
P|, is obtained only by extrapolation from the P| > region. It can be done in 
the following way. First, we rewrite the Euclidean BS equation in the form 

Mq;Pe)= [ M AB {q E -k E -P E )(t> B {k-P E ) (4.54) 

Jk E 

where <pA and <pB denotes the set of amplitude. This equation should not have 
solution at P E > because, if it exists, it is a tachyon solution. Therefore we 
instead solve the eigenvalue equation 

\(P e )MQe] Pe) = I M AB {q E - k E - P E )Mk E ; P E ) (4.55) 

J k 

for a fixed P E > 0. Then we extrapolate the eigenvalue to P E < as a function 
of P E and search for the on-shell point A(— M B ) = 1. We employ the quadratic 
function of P E for the extrapolation. We also check the results using the PCAC 
relation. 

4.6 Normalization and Decay Constants 
4.6.1 Normalization 

In order to obtain the decay constants of the pseudoscalar mesons, it is necessary 
to normalize the BS amplitude derived from the inhomogeneous BS equation. 



4.6. Normalization and Decay Constants 
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The normalization condition of the BS amplitude is given by 



1 



K f(-q 2 +)f(-q 2 ) XnimiK "~ ' dP» 

The explicit form of this equation is shown in the appendix. 
4.6.2 Calculation of the Decay Constant 

Using the normalized BS amplitude, the decay constants of the pseudoscalar 
mesons are obtained by 



R 



r B Jq 



/(-<£)/(-£) 

xTr 



2 



x R (q; p b ) Jy H - J ' JK H+I r+ (f(-ql) - f{-q 2 -M 



+E a (q ] P) + F a (q;P)}} 



(4.57) 



+ \K, 



-K 



n nmm 



(4.58) 



Kn(q;P) = -2G Dl5 lcoior j k i /( _ fe2) 1 /( _ F) tr (gC) [^ 92/2 (fe)] 

x ^ _ |)2 _ ^2 _ ^ _ w (_ q 2_ 3p2_ A .2_ Z 2 + p ^| 

Xe ff2f3 e 99293^ T (DC) 



S£(Z)iy 



- 1 9a /3 



+ j w + ^)2 _ fc 2 _ ^ _ w ^_ q 2 _ 3 p2 _ k 2 _ j2 _ p{\ | 



X 



J 93/3. 



(4.59) 
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The terms E a (q;P) and F a (q;P) represent the corrections to the Noether axi- 
alvector current due to the momentum dependencies interaction. 

For the rj and rj' mesons, the flavor structure of the BS amplitude depends 
on the relative and total momenta in general. Therefore we can not fix A" in the 
definition of the decay constant from the flavor structure of the BS amplitude. 
Instead we only consider the decay constants associated with the octet and singlet 
axialvector currents for the rj, r/' mesons, i.e., we define and calculate /§ , /q ,/§ 
and fif. 

Since the BS equation is homogeneous, the overall sign of the BS amplitude, 
and therefore the decay constants, cannot be determined. We choose the sign as 

f /* > for rj 
\ /°, > for rf . 

4.7 PC AC relation 

In our approach, there are eight axial- vector currents, JgJa = 1, • • • , 8) , which 
satisfy the PCAC relation 



d^J^(x) = 2 [m J 5 } a (x) (4.60) 

r \ J- /(d 2 )moA Q + \ a m f(d 2 ) 

[moM (x) = ^«7s ^ (4.61) 

We take the matrix element of both the sides of eq. (4.60) between a meson 
state (Pb\ and the vacuum |0) and obtain 

-r B M 2 B = 2[m S%] (4.62) 

where j% means the decay constant f w , /§ , /§ . [m ^] is defined by 



[m Q £%\ = lim % I ^ ( 1 + 1 - ) tr 



X R (q; P)m 7 5 y 



(4.63) 



p^p b J q 2\f{-ql) f{-q%) j 

Since this relation must be kept on the mass shell, we can use this as a check of 
the obtained mass by using the extrapolation. We introduce the ratio R defined 
by 



4.7. PCAC relation 
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As this ratio should become unity at P| = — M|, the value of R a (—Mg) is a good 
measure of the accuracy of the calculation, especially the extrapolation from the 
P| > region to Pf = — M\ point. 
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Chapter 5 

Results and Discussion 



In this chapter, we show the numerical results and give the disscussion. 

5.1 Parameters 

In the present approach, there are eight input parameters: quark mass m qQ for 
the up and down quarks, the current quark mass m s0 for the strange quark, 
the scale parameter of QCD running coupling constant Aqcd, the infrared cutoff 
tiF, the smoothness parameter to, the strength parameter of the KMT Go, the 
parameter of the weight function of the KMT k and the ultraviolet cutoff A uv . We 
take the ultraviolet cutoff A uv = 100 [GeV] and at such large value the physical 
observables do not depend on it. The quark masses m qQ ,m s o are chosen so 
that the renormalized masses at the momentum scale \x = 2 [GeV] become the 
m qR = 4.5[MeV] and m sR = 150[MeV], respectively. The k parameter is taken as 
k = 0.7[GeV" 2 ], which corresponds to the instanton of the average size 1/3 [fm] 
as shown in Sec. [4.1]. The other parameters Aqcd, £if, to, Gd are chosen as 
the pseudoscalar meson masses M n , M v , M v > and the pion decay constant f n are 
fitted to the experimental values. The parameters we used are A QCD = 600 [Me V], 
t = -6.89, tjp = 0.204 and G D = 75[GeV" 5 ]. 

The results of the calculation with these parameters are M n = 136 [Me V], 
M v = 515 [MeV], M v = 982 [MeV] and f n = 95[MeV]. These are in agreement 
with experimental values in less than 6% of deviation. Table 5.1 summarizes all 
the values of the parameters. 

Although Aq CD is somewhat larger than the standard value A QC d = 100 ~ 
300 [MeV], a large Aqcd is necessary in the ILA to generate the desired dynamical 
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m qR (2[GeV}) 


4.5 [Me V] 






m sR (2[GeV]) 


150 [MeV] 






Aqcd 


600 [MeV] 




M« 


136 [MeV] 


tip 


0.204 




Mr, 


5 15 [MeV] 


to 


-6.89 






982 [MeV] 


G D 


75[GcV -5 ] 




u 


95 [MeV] 


K 


0.7[GcV -2 ] 






Auv 


100 [GeV] 







Table 5.1: The values of the the parameters in the ILA approach and the obtained 
observable quantities. 

chiral symmetry breaking (D%SB). It may indicate that the chiral symmetry 
breaking is not fully explained by the ladder approximation of the one gluon 
exchange with {7a (1) breaking effect. If we take a strong KMT interaction, we 
may reduce Aqcd, but then 77 — 77' mass splitting becomes too large. As a result, 
we should choose a large Aqcd so that the D^SB is generated mainly by the 
gluon exchange interaction. Although the results of our calculation show that the 
contribution of the KMT interaction to the quark mass function and to the quark 
condensate is small (< 10%), it gives significant effects on the mass spectrum and 
the mixing angle. Then this KMT interaction can be regarded to induce strong 
{7a (1) breaking. 



5.2 Solution of the SD equation 

Let us now discuss the solutions of the SD equation. Our numerical results of 
the running quark masses are shown in Fig.5.1 and 5.2. The values of the quark 
condensates and mass function at P| = are given in the Table 5.2. 

One sees that the mass function depends weakly on Go- When Go increases 
from to 75 [GeV -5 ], B q (0) increases only by 15%. It becomes clearer by com- 
paring the effects in the SD equation of the one gluon exchange term with that 
of the KMT term. It is shown in Figs. 5.3 and 5.4. 

The dependence of the condensates of u, d quark on Go is also small as shown 
in Fig. 5.5. When G D increases from to 75 [GeV -5 ], — (qq)]( 3 changes by 13 
%. But at the same time (qq) R changes by 24%, and indicate that the effects of 
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Figure 5.1: q\ dependences of the mass function of the q quark B q with Gd = 
0,25,50,75, 100[GeV~ 5 ] 



G D [GeV 5 ] 


£ 9 (0)[GeV] 


£? s (0)[GeV] 


-(qq)T[GeV] 


-< S - S )^[GeV] 





0.524 


0.744 


0.225 


0.052 


25 


0.551 


0.752 


0.236 


0.088 


50 


0.582 


0.763 


0.248 


0.122 


75 


0.616 


0.778 


0.259 


0.146 


100 


0.655 


0.797 


0.270 


0.165 



Table 5.2: Dependences of the mass function at q\ — and the quark condensates 
on the strength of the KMT interaction. 



the KMT interaction is not so small. The dependence of the condensates of s 
quark on Gd is rather large. When Gd increases from to 75 [GeV~ 5 ], — (ss)]( 3 
changes by 64 %. 

These results clearly differ from the result of NJL model. In the NJL model, 
the effect of the KMT interaction in the dynamical quark mass are very large, 
Gd (ss) /Gs — 0.44 suggests that the contribution from the KMT interaction to 
the dynamical quark mass is 44% of that of the four quark interaction. This 
difference will be examined later. 
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0.8 



0.6 

B s [GeV] 

0.4 

0.2 



0.5 1 1.5 2 

q E 2 [GeV 2 ] 

Figure 5.2: q 2 E dependences of the mass function of the q quark B s with Go = 
0,25,50,75, 100[GeV~ 5 ] 

5.3 Solution of the BS equation of the Pseudo 
Scalar mesons 

Let us now turn to the discussion of the solutions of the BS equation. Our 
numerical results for the pion are summarized in Table 5.3. From these results, 
it is shown that the pion mass is not so sensitive to G D , which suggest that the 
increase of dynamical quark mass and of the attractive KMT interaction tend to 
cancel with each other. The change of decay constant is also small. It is 7% and 
this result is consistent with the change of the condensate. The R value defined 
in Eq.(4.64) is an indicator which shows whether the Ward identity of the axial 
vector current is satisfied on the mass shell. We consider that the deviation of 
R from the unity indicates error in the extrapolation from Euclid (p 2 E > 0) to 
Minkowski (p 2 E < 0) region. From the Table 5.3, one sees that the deviation is 
about ~ 6% and thus extrapolation error for M 2 and f n is of this order. 

The masses and decay constants of the rj and rj' are shown in Tables 5.4, 
5.5 and Fig. (5.6). Since the BS equation is homogeneous, the sign of the BS 
amplitudes, and therefore the decay constants, is undetermined. We choose the 
sign so that f 8 is positive for n, and f is positive for n'. 

The masses of r\ and n' mesons and their decay constants depend strongly on 
the Go- Especially, the n' meson mass seems to be sensitive as is expected. This 
is in contrast to the result for the pion. The ^7^(1) breaking effects are large in 
the r) and rj' sector. The overall behavior of the mass spectrum agrees well with 
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Figure 5.3: Effects of the gluon exchange interaction and the KMT interaction 
to the q quark mass function at Go — 75[GeV~ 5 ] 



G D [GeV~ a ] 




il4[MeV] 


U[MeV] 


R 





126 


84 


1.00 


25 


130 


88 


1.02 


50 


133 


92 


1.03 


75 


136 


95 


1.05 


100 


139 


99 


1.06 



Table 5.3: Dependences of the solution of the pion BS equation on the strength 
of the KMT interaction. 

the NJL result. [19, 4] 

The 77 — 77' mixing angles are defined in terms of the decay constants as 

^ = tan0„ (5.1) 

J8 

: C = tan^. (5.2) 
Jo 

The results are presented in Tables 5.4, 5.5 and Fig. (5.7). In the limit of G D — 0, 
since there is no flavor mixing, 77 is pure (uu + dd) state and 7/ is ss state. 
They are the ideally mixed states, i.e., 9 = arctan(— y/2) = —54.7°. In this limit, 
as Gd becomes larger, the mixing angle approaches to 0° i.e., the SU(3) limit. 
At that time, 77 approaches to the pure octet -^(uu + dd — 2ss) state and 77' 
approach to the pure singlet 4= (uu + dd + ss) state. At G D = 75 [GeV -5 ], we 



total mass function B q (q E 2 ) 




j 1 1 i ] 
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G D [GeV- b ] 


M„[MeV] 


/, 8 [MeV] 


/°[MeV] 


R s 


mixing angle [deg] 





134 


49 


69 


1.00 


-54.7 


25 


311 


61 


65 


1.06 


-46.8 


50 


409 


78 


53 


1.08 


-34.2 


75 


447 


95 


35 


1.08 


-20.0 


100 


448 


106 


16 


1.03 


-8.6 



Table 5.4: Dependences of the solution of the 77 meson solution of the coupled 
channel BS equation on the strength of the KMT interaction. 



G D [GeV- b ] 


M v ,[MeV] 


/£[MeV] 


/£[MeV] 


R* 


mixing angle [deg] 





785 


-91 


74 


1.55 


-54.8 


25 


846 


-71 


91 


1.69 


-35.9 


50 


867 


-54 


103 


1.71 


-34.2 


75 


982 


-28 


109 


1.73 


-26.2 


100 


1273 


-11 


82 


2.05 


-6.6 



Table 5.5: Dependences of the solution of the n' meson solution of the coupled 
channel BS equation on the strength of the KMT interaction. 
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Figure 5.4: Effects of the gluon exchange interaction and the KMT interaction 
to the s quark mass function at Ge> = 75[GeV~ 5 ] 

obtain 6 V = —20.0° and 9 rj i = —26.2°. These large changes of the mixing angle 
and the mass spectrum at Go = 75[GeV~ 5 ] indicate strong Ua(1) symmetry 
breaking effect. 

Finally we check the extrapolation ambiguities in the r\ and n' meson masses. 
The R values for the n meson are as good as the pion case. On the other hand, 
the deviation for the n' meson is as much as 73%. As the n' mass is much larger 
than the tt and r\ masses, the extrapolation from the Euclid region to the on shell 
momentum brings large uncertainty. It is shown in Fig. 5.10. The obtained mass 
from the condition R — 1 is 814[MeV] and the mass obtained from the condition 
A = 1, where A is the eigenvalue of the BS equation, is 982 [Me V]. This difference 
of the masses may show the extrapolation uncertainty. Since the dependence of 
the decay constants f®, and on the momentum P| is small, the uncertainty of 
the decay constants is smaller than that of the masses. Those are shown in Fig. 
5.10. We think that this ambiguity from the extrapolation is unavoidable in the 
current approach. 

5.4 Solution of the BS equation of the Scalar 
Mesons 

In this section, we present the results of the ILA calculation for the scalar mesons. 
Note that the parameters of the present approach are completely fixed in the 
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Figure 5.5: Dependence of the quark condensates on the strength of the KMT 
interaction. 

study of the pseudoscalar mesons. The results for the a , cr, f mesons are 
summarized in Tables 5.6, 5.7 and 5.8. The dependence of the mass spectrum on 
the strength, Go of the KMT interaction is shown in Fig. 5.11. 

First, the dependencies of the masses of a and a on G D look qualitatively- 
same as the NJL results shown in Fig. (3.1). In the NJL calculation, the parame- 
ters are chosen so as to reproduce the M w and f w at each G D . In contrast, in the 
ILA approach we change Gjj independently. However, since M n and /„. depend 
weakly on G D) the results of ILA show similar behavior as those of NJL. We note 
that the mass of a grows as Go increases, while the a mass is fairly constant. 
The behavior of the f mass differs from the that of NJL model. As noted in 
Chapter 3, the mass of the f solution in the NJL model is much higher than the 
threshold, 2M q and the momentum cutoff,A. The f mass in the NJL result is 
not so reliable. 

The a meson mass is predicted as about 670 [Me V]. This rather small a meson 
mass is interesting. In the case of the NJL model, the a meson mass is determined 
to be close to twice of the dynamical quark mass. On the other hand, in the ILA 
approach the value of the mass function at P|, B q (p 2 E = 0) is about 6 16 [Me V], 
which is comparable to the a meson mass. Although there are such differences, 
the properties of the physical observables agree in both the calculations. This 
feature is very interesting. 

The dependence of the a meson mass on the G D is small. This means that 
the effects of the ss mixing and the KMT interaction are balanced. On the other 
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Figure 5.6: Dependence of the mass spectrum of the n, r\ and r( mesons on the 
strength of the KMT interaction. 



hand, in the case of ao, the flavor content does not change and only the repulsion 
of the KMT is effective. Thus the a mass monotonically grows as Go increases. 
We obtain M ao = 942 [Me V]. This result is comparable to the experimental value 
984.8 ± 1.4[MeV]. We obtain a significant mass splitting between the a and a , 
about 272 [Me V]. We conclude that the observed a — a mass splitting can be 
explained as the ^7^(1) symmetry breaking effects. 

The obtained mass of f is larger than the mass of a by about 394MeV 
and they do not become degenerate even if we change the strength of the KMT 
interaction in our parameter range. Therefore /o(980) may not be identified with 
the qq state. To reproduce / (980) and the higher 1 = scalar meson states, it 
may be needed to treat mixings of multi-quark states. 

Next, we consider the mixing angles. We introduce the matrix elements S* 8 
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G D [GeV~ b 




M ao [MeV] 





591 


25 


694 


50 


808 


75 


942 


100 


1128 



Table 5.6: Dependences of the solution of the do BS equation on the strength of 
the KMT interaction. 



G D [GeV~ b 




M a [MeV] 


S° 


S H 


mixing angle [deg 







591 


12.9 


9.09 


-54.7 


25 


618 


13.5 


8.40 


-58.1 


50 


645 


14.4 


7.51 


-62.5 


75 


667 


16.0 


6.47 


-68.0 


100 


683 


18.1 


5.25 


-73.8 



Table 5.7: Dependences of the solution of the ameson solution of the coupled 
channel BS equation on the strength of the KMT interaction. 



G D [GeV~ b 




M /0 [MeV] 


S° 


S* 


mixing angle [deg 







1205 


-21.1 


14.9 


-54.7 


25 


1230 


-20.8 


7.4 


-70.5 


50 


1267 


-28.4 


3.1 


-83.7 


75 


1336 


-35.5 


3.8 


-83.9 


100 


1462 


-45.9 


10.3 


-77.3 



Table 5.8: Dependences of the solution of the /o meson solution of the coupled 
channel BS equation on the strength of the KMT interaction. 
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Figure 5.7: Dependence of the mixing angle of the rj and rj' meson on the strength 
of the KMT interaction 



and S° which are defined by 



, ^8 

S a = J d 4 x(0\q—q(x) (scalar meson state) 



= tr 



(5.3) 
(5.4) 



Since these S values extract the particular flavor component of (ps which is the 
main component of the BS amplitude at the origin, we employ S 8 and S° to 
determine the ratio of the octet and the singlet components. Accordingly we 
define the mixing angles of the scalar mesons as 



tan# a = 



cfo 

= 

f° ~ q/o' 



tan^ 



(5.5) 
(5.6) 



The results are summarized in Table 5.6, 5.7 and Fig 5.12. 

The mixing angle of a approaches —90°, where a becomes the purely flavor 
singlet state: -^(uu + dd + ss) . The obtained angle at G D — 75[GeV~ 5 ] is —68.0° 
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Figure 5.8: The extraporation of the eigenvalue A of the BS equation and R 
value. The solid line shows the extraporation of the A and the doted line shows 
the extraporation of the R. 

Figure 5.9: Extrapolation of the eigenvalue and the R for rj' meson with Go = 
75[GeV~ 5 ]. 

and is slightly smaller than the result of the NJL model, —77.3°. This angle 
corresponds to about 5% mixing of the strangeness component in o. 

Finally , the comparison of the scalar and pseudoscalar meson spectra is shown 
in Fig. (5.14) 

5.5 Solution of the BS equation of the Strange 
Mesons 

In this section, we calculate the masses of the strange pseudoscalar meson K and 
the strange scalar meson Kq. Here we employ a crude approximation in order 
to avoid the technical difficulty As these mesons consists of the a strange quark 
and a nonstrange quark. Then the BS equation of the ILA approach becomes 
extremely complicated. Instead of treating the asymmetric BS equation, we solve 
the symmetric BS equation for the quarks whose mass is 77.25 [MeV]. This mass 
is the average of m q and m s . The results are summarized in Tables 5.9, 5.10 and 
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Figure 5.10: Extraporation of the /°, and with G D = 75[GeV~ 



G D [GeV" & ; 




M x [MeV] 





547 


25 


546 


50 


544 


75 


538 


100 


530 



Table 5.9: Dependences of the solution of the K meson BS equation on the 
strength of the KMT interaction. 



Fig.(5.14). 

Again the K mass depends little on Go as in the case of ix . On the other 
hand, the Kq meson mass increases due to the repulsion of the KMT interaction. 

Although these results are based on the symmetric BS equations, the obtained 
K mass, 494 ~ 498 [MeV], reproduces the observed value and Go dependence 
agrees with that of the NJL. 

One sees that the slopes of the a mass and Kq mass are nearly equal and 
therefore it is clear that the Kq mass does not become smaller than the a mass 
in this approach. Since the NJL also gives the similar result, we conclude that 
the reversal of the k(700 — 900) mass and the a (980) mass can not be explained 
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Figure 5.11: Dependence of the mass spectrum of the a , cr and / mesons on the 
strength of the KMT interaction. 



by the Ua(1) symmetry breaking. Furthermore, our Kq mass is smaller than the 
next excited state -fTo(1430) b y about HOMeV. Consequently, the light 1 = 1/2 
scalar mesons may not be explained without considering mixings of the multi- 
quark states. 
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G D [GeV 5 ] 



Figure 5.12: Dependence of the mixing angle of the n and rf meson on the strength 
of the KMT interaction. 



G D [GeV~*[ 




M K *[MeV] 





946 


25 


1040 


50 


1158 


75 


1321 


100 


1542 



Table 5.10: Dependences of the solution of the Kq meson BS equation on the 
strength of the KMT interaction. 
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Figure 5.13: Dependence of the mass spectrum of the n, r], r] f , ao, o and fo mesons 
on the strength of the KMT interaction. 
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Figure 5.14: Dependence of the mass spectrum of the scalar and pseudoscalar 
meson nonet on the strength of the KMT interaction. 
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Chapter 6 
Conclusion 



We have studied the roles of the Ua(1) breaking interaction of QCD in the spec- 
trum of the light scalar nonet mesons using the extended NJL model and the 
ILA of QCD. We have first analyzed qualitative properties of the C/a(1) breaking 
interaction in the extended NJL model, and next using the ILA approach we 
have carried out quantitative study. Since these approaches are consistent with 
chiral symmetry, the scalar and pseudoscalar mesons can be regarded as chiral 
partners. We choose parameters to reproduce the masses and decay constants of 
the pseudoscalar nonet mesons and apply those to the scalar nonet mesons. 

The extended NJL model is the SU(3) NJL model supplemented with the 
Kobayashi-Maskawa-'t Hooft (KMT) interaction, which causes the C/a(1) break- 
ing. The strength of the KMT interaction is chosen so as to explain the electro- 
magnetic decays of the rj meson. To study the role of the Ua(X) breaking, we 
study the mass spectrum and the mixing angle of the scalar mesons as functions 
of the strength of the KMT interaction. As a result, we observe that the a meson 
mass changes little, which suggests that the increase of the dynamical quark mass 
and the attractive KMT interaction cancel with each other. On the other hand, 
the mass of ao increases monotonically due to the repulsive KMT interaction. 
This a — a mass splitting amounts to a 150 MeV and is slightly smaller than 
the current experimental data. A possible reason that the ao mass is not large 
enough is the large unphysical qq decay width which is an artifact of this model. 
We have also found that the strangeness content in the a meson is about 15%. 

The ILA of QCD is an approximation that is consistent with chiral symmetry 
and consists of the rainbow approximation of the Schwinger-Dyson equation and 
the ladder approximation of the Bethe-Salpeter equation. Using this approach, 
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we analyze the scalar meson spectrum quantitatively. As we do in the extended 
NJL model, we choose the parameters to reproduce physical observables of the 
pseudoscalar mesons. 

We have obtained that the a meson mass comes around 600 ~ 700 [Me V], 
which is much lower than the value expected from the effective quark mass. This 
suggests that the a meson is special as the chiral partner of the pion. We obtain 
the a — do mass splitting of about 272 [MeV] and the strangeness content in the a 
meson is about 5%. It will be interesting to check this ss mixing experimentally 
by using, for instance, the radiative J/ip decays in near future. 

The obtained mass of f is larger than the mass of a by about 394MeV 
and they do not become degenerate even if we change the strength of the KMT 
interaction in our parameter range. Therefore /o(980) may not be identified with 
the qq state. To reproduce /o(980) and the higher 1 = scalar meson states, it 
may be needed to treat mixings of multi-quark states. 

The masses of the strange mesons have been calculated in an approximation. 
We have found that the Kq mass is larger than the a mass by about 230 [MeV]. 
The masses of a and Kq mesons increase by almost the same rate as Gd is 
increased. Therefore it seems that the reversal of the k(700 — 900) mass and the 
ao(980) mass can not be explained by the Ua(1) symmetry breaking. 

Furthermore, our Kq mass is smaller than the next excited state Kq(1A30) 
by about 379MeV. Consequently, the light 1 = 1/2 scalar mesons may not be 
explained without considering mixings of multi-quark states. 

In conclusion, we have shown that the spectra of the pseudoscalar meson 
nonet and the a and a masses can be reproduced by the extended Nambu-Jona- 
Lasinio model and improved ladder approximation of QCD with C/a(1) breaking 
interaction. We have found that the observed spectrum is consistent with the 
strong Ua(X) breaking interaction. Furthermore we have shown that / (980) and 
k(700 — 900) can not be reproduced in our picture. It may be needed to treat 
mixings of multi-quark states. 

For the more detailed analysis, the approach in which the decay width of 
the mesons may be needed. Such calculations will require two-point correlation 
functions of mesons in the Minkowski kinematics. This is beyond the scope of 
the present approach. 

As these results have clarified the realization of chiral symmetry in the scalar 
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mesons, it will now be interesting to confirm the results experimentally and also 
to study behavior of the scalar mesons at finite temperature and/or density where 
we expect that chiral symmetry is going toward restoration. 
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Appendix . Conclusion 



Appendix A 
BS Equation 



Here we show the explicit form of the BS equation. At first, we define the 
regularized amputated BS amplitude x R (q; P) by 



X R (q; P) = Sr\q + ^)x R (q; P)S§T\q - f ) 



(A.l) 



By introducing the 4>a which denotes the set of the amplitudes of the P) , 
the BS equation can be written as 



(A.2) 



0(g;P)= K AB (q,k;P)Mk;P)- 

J k 

We solve the BS equation by dividing to these two part. 
The explicit form of Eq. (A.l) is given by 



0s = -r 



1 

A 



0Q 



1 

A 



1 

A 



{-q E 2 + - B + B_} j>s 

-q E 2 (B + -B_) + ^(B + + B_)}$ P 
'£(B + + B_) - P E q E (B + - B_)} Q 
+ {q E 2 PE 2 -(P E qE) 2 }0T 

{B + -B_}4>s 
+ {-q E 2 -Zt-B + B_}^ P 

+ {-2P E q E }4> Q 

+ {(P E q E )(B + + B_)-Zf(B + -B_)}4>T 
\PEqE)<pp 

(q E 2 + ^)-B + B_}4> Q 
-q E 2 (B + + B_) + ^p(£?+ - £_)} T 



(A.3) 



(A.4) 



(A.5) 
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4>s 

\{B+-B_)}4> P 



+ {B + + 
+ {q E 2 -^-B + B^}^ T 



(A.6) 



where 



P 



qE- 



P 

q ~ 2 



A = (q E l + Bl)(q E l+B 2 _). 



(A.7) 
(A.8) 
(A.9) 



K 



ss 



. pp 



K, 



(E) 
QQ 



C F - 

Cp^-r 

c F - 

C F - 

Cp^'-r 

c F - 



{qE - k E ) 2 
1 



(A.10) 



+ 



k n 



(qE . kE?A ' (qE - kE) M-^)-^-^mv 



The explicit form of the £qe part of the interaction kernel K in the Eq. (Al) 
is given by 

7 2 2 
?\° 

> 2 

fx* 

7 2 2 
V 2 2 

> 2 

7 2 2 {q E -k E y 



kp 

t - kp) 2 A (q E 



1 2 (k P 

We - k E ) 2 ~ {q E -k E y {qE - kE) ■ Pe [a / 

\k\ - q\) ("%) + 2fe - fc B ) • P E ^-2(q E -k E )- k E 



(A.12) 
(A.13) 
(A. 14) 



where 



A = q 2 E P 2 - (P E ■ q E ) 2 

k q = P 2 E {q E ■ k E ) - (P E ■ q E ){P E ■ k E ) 

kp = q E {PE ■ k E ) - (q E ■ P E ){q E ■ k E ). 



The £rmt part is given by 



4? = -2AG D e^'e^'tr[S F (p)} h , h 



ghf'fh'g' 



(A.16) 
(A.17) 
(A.18) 



(A.19) 
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Appendix B 



Normalization Condition of the 
BS Amplitude 



Here we show the normalization condition of the BS amplitude explicitly. 
To reduce the expressions we use abbreviation 



B+- 


= B((q E + P E /2) 2 ) 


(B.l) 


B_ -- 


= B((q E -P E /2f) 


(B.2) 


P+E 


= (q E + P E /2) 2 


(B.3) 


P~E 


= (q E - P E /2f 


(B.4) 


$ = 




(B.5) 


$ = 


x R (q;P). 


(B.6) 



The normalization condition of the BS amplitude is given by 



I 



N c f 

1 



2 V f(q E l)f(qE 2 - 



X | -ti[$ f> E $fe] - -ti[$ f> E $f> E ] - -5_tr[$f B $] - \P E q E + J B' + ti[^ E ] 

+ \ (p E q E + ^fj B' + ti[$$f> E ] + (p E q E + ^ B' + B_tr[$<£>] 

-^tT[ty E $f> E ] - hr[$f E <S>f E ] + ^£ + tr[l»$f £ ] ~ [~P E q E + ^ B'J,t[He$\ 

-\ [~P E q E + ^fj B'M$P E <£>] + (-P E q E + £+#_tr[$$]} 
= 2P E (B.7) 
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tr[l»^$f E ] 

ti[$ f> E $f> E ] 
tr[$f £ $] 

tr[$$f B ] 

tr[$^$] 

tr[$$^] 

tr[$$] 
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4 (-M-q; p )Mq; p)Peqe - M-q; p )Mq; p^Peqe) 2 - Phi) 
-M-q; p)<Ps(t, P)((PEq E ) 2 - pWe) + M-q; p)Mq; P)q\P E q E 
+M-q; P)Mq; P)(P E q E ) 2 - M~q; p)Mq; P)q\Pl 

-M-q; P)0 Q (g; P)P 2 E P E qE) (B.8) 

4 (-M-q; p )Mq; P)P E q E + M-q; p)Mq; P)((P E q E ) 2 - Phi) 
+M-q; p)Mq; P)((P E q E ) 2 - PhD + M-q; p)Mq; p^Peqe 
+M-q; P)Mq; P)q\Pl - M-q; p)Mq; ^)(^) 2 

-</> Q (-q; P)<P Q (q; P)P 2 E P E q E ) (B.9) 

4 (-M-q; p)Mi; p)p 2 e - M-q; p)Mq; p)H{(Peqe) 2 - PhD 

+M-q; p)Mq; P)(P E q E ) 2 + M-q; P)Mq; P)P E P E q E 

-M-q; P)Mq; P)P 2 E P E q E - M-q; p )Mq; p)H) ( B -io) 
-4 (-M-q; P)Mq; P)P E q E - M-q; P)Mq; P)P E q E 

-M-q; P)Mq; p)p 2 e + M~q; P)Mq; p)p 2 e 

-M-q; P)Mq; P){(P E q E ) 2 - PhD 

-M-q; P)Mq; P){(P E q E ) 2 - Ph 2 E }) (b.ii) 
-4 (-M-q; P)Mq; P)P E q E - M-q; P)Mq; P)P E q E 

-M-q; P)Mq; p)p e + M-q; p )Mq; p)p 2 e 

+M-q; P)Mq; P){(P E q E ) 2 - Ph 2 E } 

+M-q; P)Mq; P){(P E q E ) 2 - PhD) (B.12) 
-4 (-M-q; P)Mq; P)q% - M-q; p)Mq; p)q% 

-M-q; P)Mq; P)P E q E + M-q; p)Mq; P)P E q E 

+M-q; P)Mq; P){(P E q E ) 2 - PhD 

-M-q; P)Mq; P){(P E q E ) 2 - PhD) (b.13) 

-4 (-M-q; P)Mq; p)<& - M-q; P)Mq; p)<& 
-M-q; P)Mq; P)P E q E + M-q; p)Mq; p)Peqe 
-M-q; P)Mq; P){{P E q E ) 2 - PhD 
+M-q; P)Mq; P){(P E q E ) 2 - PhD) (B.i4) 

4 (M-q; P)Mq; p) - M-q; p)Mq; p)q% - M-q; P)Mq; p)Peqe 
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+0 Q (-g; P)M<i; P)P E q E + </> Q (-q] P)M^ P)P E 

-M-q; p)Mq; ^H(^) 2 - PhD) ■ 

(B.15) 



1500 r- 

^ 1400 - 
> 1300 - 

CD 

^ 1200 - 

C/> 1100 - 
CD 

<g 1000 - 
[j? 900 - 
800 - 
700 - 
600 - 
500 1 
400 - 
300 - 
200 - 
100 - 



25 50 75 100 

G D [GeV 5 ] 



300 




G D [GeV" 5 ] 



1300 




co 700 - 

CO 

gj 600 - 
500 - 
400 - 
300 - 
200 - 
100 - 



K 



25 50 75 100 

G D [GeV 5 ] 



1300 
1200 
1100 
1000 
900 
800 



m 700 

CO 

600 



500 
400 
300 
200 
100 






G D [GeV J ] 



